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Preface 


The  purpose  of  this  effort  was  to  determine  if  vortex 
paneling  methods  can  compute  the  vortex  bursting  that  occurs 
over  a  unit  aspect  ratio  delta  wing  at  high  angles  of 
attack.  The  leading-edge  vortex  flows  dominate  the  flow 
field  over  such  wings.  The  vortex  breakdown,  or  bursting 
that  occurs  over  the  wing  affects  the  lift  and  drag,  as  well 
as  the  moment-balance  of  the  aircraft.  It  is  therefore 
desirable  to  compute  such  effects  if  possible. 

Multiple  calculations  were  performed  on  both  the  ELXSI 
at  the  Air  Force  Institute  of  Technology,  and  the  Cray-YMP 
at  Ohio  State  University.  The  unreliability  of  the  ELXSI 
provided  great  anxiety  as  to  whether  or  not  this  thesis 
would  be  completed.  However,  the  Cray-YMP  worked  quite  well 
and  was  an  invaluable,  albeit  limited,  resource. 

I  would  like  to  take  this  time  to  give  appreciation  to 
others.  I  would  like  to  thank  my  thesis  advisor  Major  Curt 
Mracek  for  his  invaluable  time,  support,  and  assistance.  I 
would  also  like  to  thank  the  other  committee  members.  Cap¬ 
tain  Phil  Beran,  and  Lieutenant  Colonel  DeJongh  for  their 
valuable  inputs.  I  would  also  like  to  give  thanks  to  all  my 
classmates  for  their  help  and  friendship.  Finally,  and  most 
importantly,  I  would  like  to  give  thanks  to  my  wife,  Lora, 
for  her  support  and  understanding  throughout  the  past  one 
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and  a  half  years. 
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Abstract 


A  unit  aspect  ratio  delta  wing  is  analyzed  to  determine 
if  a  unsteady  vortex  panel  method  can  calculate  the  vortex 
bursting  of  the  vortex  developed  along  the  leading-edge.  A 
two-dimensional  vortex  panel  method  is  used  to  investigate 
wake  splitting  schemes  of  the  vortex  cores  that  comprise  the 
vortex  sheet,  or  wake.  The  worthier  of  the  splitting  meth¬ 
ods  is  later  implemented  into  the  three-dimensional  wake. 

The  weighting  values  of  the  Kutta  equations  are  investigated 
for  the  leading-edges,  and  are  determined  to  be  problem  spe¬ 
cific.  An  alternative  approach  to  the  management  of  the 
wake  data  is  presented  and  implemented.  The  delta  wing  is 
then  analyzed  for  20.5  and  30  degrees  angle  of  attack.  The 
panel  method  produced  a  perturbation  and  enlargement  in  the 
wake  at  the  two-thirds  chord  position  over  the  winq  at  30 
degrees  angle  of  attack,  indicative  of  vortex  breakdown. 
However,  a  qualii'ative  comparison  with  the  results  of  the 
20.5  degree  angle  of  attack  calculation  negates  this  conclu¬ 
sion.  Vortex  bursting  was  not  evident  over  the  wing  at  Ju 
degrees  angle  of  attack. 


chapter  I 

Introduction 

1.1  Problem  Statement 

To  determine  if  vortex  paneling  methods  can  compute  the 
vortex  breakdown  that  occurs  over  a  delta  wing  at  high 
angles  of  attack. 

1.2  Motivation  for  Research 

The  design  of  modern  aircraft  combines  both  analytical 
and  experimental  testing.  As  more  powerful  computers 
emerge,  the  usefulness  and  applicability  of  computer  codes 
to  predict  the  fluid  characteristics  over  airfoils  has 
increased.  As  aircraft  operate  in  the  high  angle  of  attack 
regimes,  the  flow  field  becomes  quite  complex  and  is  domi¬ 
nated  by  vortex-flow  effects.  The  aerodynamic  challenges 
this  brings  are  expected  to  increase,  as  designers  try  to 
solicit  the  advantages  of  this  flow.  The  ability  to  compute 
this  complex  airflow  would  facilitate  the  design  process  and 
reduce  the  costs  associated  with  development. 

1.3  Background 

A  grasp  of  what  is  vortex  breakdown  and  how  will  it  be 
determined  numerically  needs  to  be  addressed.  Within  the 
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context  of  this  effort,  the  terms  "vortex  breakdown"  ana 


"vortex  bursting"  are  used  synonymously  and  are  restricted 

in  reference  to  the  fluid  flow  over  delta  wings.  Ekaterina 

ris  and  Schiff  describe  vortex  breakdown: 

An  increase  in  the  angle  of  attack  strengthens  the 
vortices  (generated  along  the  leading  edges)  until 
eventually  a  sudden  change  occurs  in  the  nature  of 
the  cores.  This  sudden  change  is  know  as  vortex 
breakdown.  Vortex  breakdown  involves  a  transition 
of  the  vortex  core  from  jet-like  to  a  wake-like 
flow.  (2:1) 

Hitzel  defines  vortex  breakdown  as  "...  a  rapid  change  of 
the  vortex  structure  and  a  corresponding  diminishing  of  the 
induced  suction  forces"  (3:73).  O’Neil,  Barnett,  and  Louie 
describe  the  breakdown  effect  with  respect  to  the  leading- 
edge  vortex,  "...  the  well-defined  leading-edge  vortex 
degenerated  rapidly  into  a  substantially  larger,  more 
diffuse  vortical  flow  region  with  relatively  mild  gradients 
(7:220) . 

Several  attempts  have  been  made  to  compute  the  leading 
edge  vortex  breakdown  that  occurs  over  delta  wings  at  high 
angle  of  attack.  Hitzel  investigated  this  phenomena  using 
Eul er-methods .  He  found  that  the  Eul er-methods  simulated 
the  leading-edge  vortices  despite  the  neglect  of  viscous 
effects.  He  states. 

Often  vortex-breakdown  was  supposed  to  be  triggered 
by  viscous  effects.  However  as  experiments  showed, 
breakdown  is  almost  independent  of  viscous 
ef fects ...  Breakdown  is  triggered  by  adverse 
pressure-gradients  which  decelerate  the  axial  flow. 

By  the  demands  of  continuity  and  the  conservation  of 
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momentum  the  vortex-core  is  forced  to  widen  conside¬ 
rably... and  the  suction-peaks  producing  the  nonlin¬ 
ear  lift  will  decline. 

In  his  calculations,  Hitzel  was  able  to  compute  the  leading- 
edge  vortex  breakdown  that  provided  "good  agreement"  with 
experiments.  However,  his  solutions  "exhibit  very  high 
total  pressure  losses  in  the  vortex-core",  and  he  attributes 
the.se  losses  to  numerical  errors  introduced  by  the  descreti- 
zation  of  the  mesh.  Furthermore,  there  were  two  vortex 
diminishing  effects.  One  of  these  causes  was  the  decreasing 
geometrical  resolution  of  the  vortex  wake.  As  the  vortices 
leave  the  wing,  they  flow  into  the  coarser  outer  part  of  the 
computational  domain.  The  other  vortex  diminishing  effect 
was  caused  by  the  downstream  boundary  conditions.  Hitzel 
states,  "Any  boundary  condition  impose  some  constant  free¬ 
flow  conditions,  (therefore)  the  vort  t  almost  must  cease  to 
exist  at  the  boundary."  (3:73-83) 

O'Neil,  Barnett,  and  Louie  also  investigated  vortex 
breakdown  over  delta  wings  using  Eul et -methods .  They  echo 
Hitzel 's  remarks  that  vortex  breakdown  is  governed  primarily 
by  inviscid  factors.  They  found  in  their  computations  that 

the  well-defined  leading-edge  vortex  degenerated 
rapidly  into  a  substantially  larger,  more  diffuse 
vortical  flow  region  with  relatively  mild  gra¬ 
dients... The  occurrences  of  such  breakdown- 1 ike 
events  in  the  Euler  solutions  were  found  to  have 
definite  trends,  and  were  clearly  not  a  result  of 
arbitrary  numerical  coincidence. 
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They  developed  a  method  to  identify  the  location  of  vortex 
breakdown.  At  the  onset  of  vortex  breakdown,  the  locus  of 
maximum  swirl  angle,  falls  in  toward  the  vortex  axis.  This 
identifies  the  point  of  vortex  breakdown,  or  shear-layer 
collapse.  However,  they  too  experienced  the  numerical  dis¬ 
sipation  induced  by  the  artificial  viscosity,  which  is  com¬ 
mon  to  all  Euler  codes.  (7:218-226) 

1.4  Scope  of  Development 

The  effort  used  the  vortex  panel  method  developed  by 
Mracek.  (6:11-100)  Panel  methods  are  used  to  describe 
inviscid,  potential -fl ow  aerodynamics.  The  method  is  a 
hybrid  of  vortex  paneling  methods,  which  describe  the  sur¬ 
face,  and  vortex  lattice  methods,  which  are  used  to  model 
the  wake.  The  advantages  over  Euler-methods  are  numerous. 
Hitzel  states, 

"In  Eul er-cal cul ations  usually  the  start  conditions 
already  introduce  vorticity...  A  proper  natural 
evaluation  of  the  vorticity  however  should  start 
from  zero  evolving  to  the  f 1 ight-vel ocity  to  simu¬ 
late  the  real  time-dependent  build-up  of  the  flow." 
(3:75) 

This  was  performed  within  the  panel  method  employed  by 
starting  the  fluid  flow  impulsively.  This  method  has  a 
continuously  developing  wake  that  xs  not  bound  by  a  computa¬ 
tional  domain,  or  mesh.  This  eliminates  the  vortex  dimi¬ 
nishing  effects  and  the  initial  induced  vorticity  over  the 


4 


surface  as  found  in  Euler-codes.  This  method  also  allows 
for  the  ability  to  apply  and  move  control  surfaces  on  the 
airfoil,  and  account  for  roll-,  pitch-,  and  yaw-rates. 

A  two-dimensional  vortex  panel  method  was  used  to 
investigate  two  separate  splitting  algorithms  of  the  wake. 
The  wake  will  be  modeled  as  a  vortex  sheet  composed  of 
discrete  vortex  cores.  As  the  wake  cores  are  convected 
downstream,  the  distance  between  adjacent  cores  i  ease. 

In  order  to  sustain  a  continuous  distribution  of  wake  cores, 
additional  cores  are  added  between  those  that  have  moved 
sufficiently  apart.  The  strengths  of  the  added  cores  are 
determined  from  the  conservation  of  circulation  within  the 
wake.  This  provided  a  splitting  model  to  be  used  in  the 
three-dimensional  vortex  panel  method.  Similar  to  the  two- 
dimensional  modeling  of  the  wake,  the  three-dimensional 
representation  of  the  vortex  sheet  used  discrete  vortex 
cores,  although,  they  are  of  finite  length  due  to  the  added 
dimension.  An  investigation  into  weighting  values  on  the 
Kutta  condition  for  the  leading  edge  was  accomplished  to 
create  the  correct  flow  condition  off  the  leading  edge.  The 
wake  data  was  organized  into  a  finite  element  methodology. 
This  reduced  the  amount  of  computer  memory  and  computational 
time  required  for  a  given  calculation.  An  analysis  is  then 
presented  on  a  unit  aspect  ratio  delta  wing,  comparing  the 
results  at  20.5  and  30  degrees  angle  of  attack. 
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Chapter  II 


Two-dimensional  Model 

2.1  Introduction 

A  two-dimensional  vortex  panel  method  was  used  to 
investigate  two  different  wake  splitting  algorithms.  The 
superior  of  the  two  models  was  then  extended  into  the  three- 
dimensional  vortex  panel  method.  A  panel  method  describes 
potential -f 1 ow  aerodynamics.  The  surface  of  the  airfoil  is 
modeled  using  segments  on  which  the  surface  vorticity  is 
linearly  varying.  The  wake  was  modeled  as  a  vortex  sheet 
composed  of  vortex  cores.  As  the  cores  are  convected  down¬ 
stream,  they  disperse  and  move  further  apart.  In  order  to 
model  a  continuous  wake,  additional  cores  are  added  between 
those  cores  that  have  moved  sufficiently  apart.  This  was 
termed  wake  splitting,  as  the  strength  of  the  new  core  is 
dependent  on  the  two  adjacent  cores. 

2.2  Two-Dimensional  Vortex  Panel  Method 

A  two-dimensional  vortex  panel  method  was  used  to 
investigate  the  wake  splitting  because  of  its  simplicity  and 
corresponding  relationship  to  the  three-dimensional  vortex 
panel  method.  In  both  cases,  a  potential  flow  model  is  used 
to  describe  the  flow.  The  velocity  is  described  as: 
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V 


(2-1) 


-  V<t) 

The  continuity  equation  for  incompressible  flow  is: 

VV  -  O  (2-2) 

Subsituting  in  for  the  velocity  potential  yields  Laplaces 
equation . 


V^(t)  -  0  (2-3) 

The  vortex  distribution  was  choosen  such  that  Laplaces 
equation  was  always  satisified,  therefore,  only  the  boundary 
conditions  needed  to  be  investigated. 

For  the  two-dimensional  model,  the  surface  of  the 
airfoil  was  modeled  using  line  segments,  or  panels,  on  which 
the  surface  vorticity  is  linearly  varying.  Time  was 
partitioned  into  time  steps  and  equated  to  the  average 
length  of  the  panels  representing  the  surface.  Both  the 
chord  and  freestream  velocity  was  nondimensionalized.  A 
wake  was  generated  off  the  trailing  edge  and  was  modeled  as 
a  vortex  sheet  composed  of  discrete  vortex  cores.  At  each 
new  time  step  a  new  trailing  edge  core  was  generated.  The 
wake  cores  were  then  convected  at  the  local  particle 
velocity  to  preserve  a  pressure  free  wake  position. 
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The  total  velocity  of  the  potential  flow  around  a  body 
can  be  written  in  terms  of  the  free-stream  velocity  and  the 
disturbance  velocity. 

V  -  (2-4) 

The  disturbance  velocity  is  the  change  in  velocity  due  to 
the  presence  of  the  physical  body  and  the  velocity  induced 
by  the  wake. 


*  ^body  ^wake 

The  no  penetration  condition  on  the  surface  was  used  to 
solve  for  the  vorticity  on  the  airfoil. 


V  • ”  0  (on  the  surface)  (2-6) 

where  ft  is  the  unit  normal  vector  to  the  surface. 

The  velocity  induced  at  any  location  in  space  by  a 
point  vortex  is  given  by:  (4:381) 


dV 


dr 

2itr 


(2-7) 


8 


Figure  1  presents  a  panel  where  the  vorticity  is  linearly 
varying . 


Figure  1.  Two-Dimensional  Vortex  Panel 

Where  the  linearly  varying  vortex  distribution  is  provided 
by: 


nm  -  G,(i-0 . 0^1 

Where  Gi  and  G2  are  the  vortex  strengths  at  the  endpoints. 
It  is  known  that 


dr  -  n(?)d?  (2-9) 

Therefore,  the  total  velocity  induced  by  a  vortex  panel  is 
determined  by  integrating  along  the  length  of  the  panel: 
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V. 


r 


ncvd^ 


2itr 


-5, 


Using  a  right-handed  coordinate  system  yields: 


Qq  -  -sinpf  +  cos(i; 

r  - 

where , 


sin  p 


y 

r 


cos  (3 


r 


Thus , 


V, 


2il 


(x-^)2  y 


f 


i 


-L 


X-? 


2ii  L(x-?)  = 


y  ci^ 


By  normalizing  the  length  of  the  panel  to  one,  the 
strength  along  the  panel  becomes: 

n(^)  -  G,(i-5)  -K  G3^ 


(2-10) 


(2-11) 


(2-12) 


(2-13) 


vortex 


(2-14) 
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Substituting  for 


ryf-(x-^)n 

L(x-^)^  +  y^J‘ 


(2-15) 


G.  p(i-^)[yr-(x-^);]^ 

2nJo  (x-^)  +  y^ 

^  2i  r‘ilizii££zJ2idf 

2it  Jo  (x  -  +  y^ 


(2-16) 


This  supplies  four  integral  equations  that  need  to  be 
solved: 


px  Cj-oGji  ^py  Cj-oGiy 


.^  =  0^2^  ■*■  ^pylG,  =  oG2y 


(2-17) 


where , 
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V 


V 


V 


px I Gj-O 


py I Gj-o 


px I c , «0 


V. 


L  f'_iizi2y_d? 


2nJo  (x  - 

»(l-^)(^-x) 

0  (x-^)2  +  y2 

y? 


±f 

2nJo 


d| 


-f‘— 

2iiJo  fx- 


(x-?)^*y2 


d? 


±r  ds 

2njo  (x-^)2  +  y^ 


(2-18) 


Solving  these  integrals  yields: 


Vp. 

Igj-O  “ 

1 

2lt 

Vp/ 

1 

O 

1 

O 

1 

2n 

Vpx 

1 

o 

1 

o 

1 

2n 

Vpy 

Igj-O  “ 

1 

(1  -x)A0  +  |lni? 


(1-x) 


In  y  A0  -  1 


X  A0  +  — In  ^ 
2 


1  -  yA0  -t-  -xln/? 


(2-19) 


where , 


A0  *  tan’’  —  I  -  tan 


•M 


In/?  «  In 


(X-  I)^-^y^ 

x^  +  y^ 


(2-20) 
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There  are  two  singularities:  x=0,  y=0,  and  x=l,  y=0.  Other 
regions  that  deserve  careful  consideration  are  when  y=0,  and 
x<0,  0<x<l,  and  x>l.  The  results  of  these  limits  are: 


x<0  y-0  A0-O 

x>l  y-0  A0-O  (2-21) 

0<x<l  y-»0*  A0-II 

0<x<l  y-»0'  A0--II 

A  simple  coordinate  transform  was  used  when  the  panel  did 
not  lie  on  the  axis.  Where, 


considered  as: 


cosy 

sin  ■< 

-  sin  y 

COS"' 

iduced  on  a 

point 

Vp 

-  V,i 

+  V. 

-  [V, 

V,]{ 

-  [Vx 

v,][ 

'I' 


(2-22) 


(2-23) 


cosy  siny 


siny  cosyJ\J 


I' 


where , 
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X 


(2-24) 


V,  -  V. 
V\,  -  V, 


V 


V 


pxlc,»0^2 


py  i  c ,-o^  2 


The  vector  dot  product  with  the  outward  normal,  required  for 
the  no-penetration  condition,  was  then  performed  on  this 
component  of  the  velocity. 

The  surface  of  the  airfoil  consisted  of  numerous  panels 
joined  together  at  their  endpoints  as  illustrated  in  Figure 
2.  The  endpoints  were  termed  nodes.  For  a  closed  surface, 
there  are  the  same  number  of  nodes  as  there  are  panels.  The 
surface  vorticity  was  constrained  to  be  continuous,  i.e.  the 
Gi  of  a  one  panel  was  equal  to  the  Gr  of  the  preceding 
panel.  Therefore,  each  node  has  a  single  vortex  strength 
(G). 


Figure  2.  Vortex  Panel  Modeling  of  a  NACA0012  Airfoil 
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The  no  penetration  condition  was  only  solved  at  the 
midpoints  of  each  panel,  or  control  points.  Each  panel  has 
an  influence  at  every  control  point.  When  the  vector  dot 
product  is  performed  with  the  outward  normal,  this  provides 
a  number  of  equations  equal  to  the  niimber  of  panels 
describing  the  surface.  This  results  in  an  influence  matrix 
[A].  The  only  unknowns  are  the  vorticity  at  each  node. 

Hence , 


[A]<G>  - 

where  is  the  perturbation  in  velocity  that  the  body 

produces  at  each  control  point. 

The  wake  consists  of  vortex  cores  that  have  been 
generated  at  the  trailing  edge,  then  convected  downstream  at 
the  local  particle  velocity.  The  strengths  of  all  the  wake 
cores  are  known  except  for  the  trailing  edge  core.  Its 
strength  was  determined  from: 

r..,.  -  0  (2-26) 


or , 


r  +  r 


wing 


wake 


0 


(2-27) 
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where , 


wake 

Tc 


r 


wing 


n-fids 


Ig,(,  =  [i]{G} 


(2-28) 


distance  between  control  points 
^circulation  of  the  remaining  wake  cores 
circulation  strength  of  the  trailing  edge  core 

-  r  -  r 

wake  wing 

-r.ake-[n<G> 


The  [/]  is  a  row  matrix  containing  the  distance  between  the 
control  points  of  adjacent  panels  and  {G)  is  a  column  matrix 
comprising  the  vorticity  strength  at  each  node. 

The  wake  also  influences  the  velocity  at  the  control 
points.  The  vortex  sheet  was  modeled  as  a  set  of  discrete 
cores.  The  velocity  induced  at  any  point  in  space  by  a 
vortex  core  is  determined  from:  (4:381) 


V 


2iir 


(2-29) 
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where  r  is  the  distance  from  the  vortex  core  to  the  point  in 
question . 

The  trailing  edge  core  also  induce  velocities  at  the 
control  points,  given  by: 


V 


Cl 


2rCTi 


(2-30) 


where , 

Fj  =  distance  from  trailing  edge  core  to  a  control  point 

The  vector  dot  product  with  the  outward  normal  vector  at 
each  control  point  can  be  performed.  The  result  is  thf 
disturbance  the  trailing  edge  core  has  on  the  velocity  at 
the  control  points. 

{c}r,  -  [V,]  (2- 

The  matrix  {c}  is  comprised  of  the  distances  from  the 
trailing  edge  core  to  each  control  point.  The  no 
penetration  condition  is  therefore, 

[A]{G}.{c}r,*  V„-rt  -  0  <2- 
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Substituting  for  the  trailing  edge  core  strength, 


or , 


[A*]<G>  -  [b*] 


(2-34) 


where , 


[A’]  =  [A]-<c>[Z]  (2-35) 

[b’l  - 

The  only  unknowns  are  the  strength  of  the  vorticity  at 
the  nodes,  however,  imposing  the  Kutta  condition  requires 
that  the  strength  of  the  vorticity  at  the  trailing  edge  is 
zero,  or  Gi=0.  Therefore,  the  number  of  equations  now 
exceed  the  number  of  unknowns.  The  method  of  solution  was  a 
reduction  in  the  square  of  the  errors.  This  is  given  by: 

[A’]^[A’]<G>  -  [A’]^[b’]  (2-36) 

This  matrix  equation  can  be  solved  using  a  number  of 
different  numerical  methods  techniques  available.  Once  the 
strengths  of  the  vorticity  are  known  at  the  nodes,  the 
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strength  of  the  trailing  edge  core  is  calculated.  The  wake 
cores  are  then  convected  at  the  local  particle  velocity  to 
produce  a  pressure  free  surface. 

2.3  Wake  Splitting  Algorithms 

The  wake  cores  are  convected  at  the  local  particle 
velocity.  The  velocity  induced  on  a  core  by  itself  is  infi¬ 
nite.  Since  this  is  unrealistic,  there  was  a  small  area 
around  each  wake  core  for  which  that  core  had  no  influence. 

Once  the  wake  cores  move  downstream,  large  distances 
are  created  between  adjacent  cores.  In  an  attempt  to  model 
the  diffusion  of  the  vorticity  and  to  sustain  a  continuous 
wake,  wake  cores  will  be  added  between  those  that  have  moved 
sufficiently  far  apart.  This  is  illustrated  in  Figure  3. 

Two  different  wake  splitting  algorithms  were  investi¬ 
gated.  In  each  case,  cores  are  divided  into  smaller  ones  as 
the  wake  forms.  The  first  method  of  wake  splitting  was 
similar  to  that  investigated  by  Mook,  et  al  (4:4).  Once  the 
distance  between  two  adjacent  cores  reached  a  prescribed 
distance,  \/2dt,  a  core  would  be  added  and  the  two  bounding 
cores  reduced  in  strength.  This  distance  was  choosen  in 
order  to  sustain  approximately  the  same  distance  between 
every  core.  The  strength  of  the  added  core  would  be  given 
by: 
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(2-37) 


r 


added 


n*  1 


and  then 


(2-38) 


If  additional  cores  were  required  to  be  added  on  both 
sides  of  a  wake  core,  the  strength  of  the  added  cores  would 
be : 


where , 


Ir  Ip 

^  added  ]  **  3  ^  4  " 


-r 

added2  ^  ^ 


(2-39) 


(2-40) 
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In  this  way,  the  total  circulation  of  the  wake  remained 
constant  throughout  the  splitting  process.  Figure  3 
illustrates  this  procedure.  A  row  of  four  nodes  are 
represented  by  the  darkened  circles  at  an  arbitrary  time,  x. 
At  x+6t,  they  are  convected  at  the  local  particle  velocity 
to  their  new  respective  positions.  Enough  distance  has  been 
created  between  the  last  three  nodes  to  require  splitting. 
This  necessitates  the  addition  of  two  wake  nodes  represented 
by  the  open  circles.  The  strengths  of  the  wake  nodes  are 
then  determined  by  the  equations  2-37  through  2-40. 


12  3  4 

r  «[- - 

V  ** 

>■  ' 

r  +  jr  - o — - O-  -  — 

I  2  3  4  5  6 


Figure  3.  Basic  Two-Dimensional  Wake  Splitting 


The  second  algorithm  was  an  extension  of  the  first. 

The  wake  cores  would  be  convected  at  the  local  particle 
velocity  to  establish  their  new  respective  positions.  If 
the  distance  between  any  adjacent  cores  was  sufficient  to 
warrant  splitting,  a  new  core  would  be  added  between  them  at 
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their  original  positions  before  convection.  The  strengths 
of  which  are  given  by  equations  2-37  through  2-40.  The 
entire  wake  would  then  again  be  convected  and  the  distances 
between  adjacent  cores  rechecked.  This  would  continue  until 
all  the  distances  between  adjacent  cores  were  within  the 
specified  limit.  This  method  is  illustrated  in  Figure  4. 


Pass  1 

T  *■  ir 


Pass  2 

T  *  6t 


- Qc— - c^: - ^ 

.  \  s  \  > 


-• - 0 - ^ 


Figure  4.  Enhanced  Two-Dimensional  Wake  Splitting 


2.4  Results 

The  airfoil  investigated  was  a  NACA0012  airfoil.  The 
following  equation  was  used  to  determine  the  location  of  the 
nodes.  (1:113) 
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A  unit  aspect  ratio  delta  wing  is  analyzed  to  determine  if  a  unsteady 
vortex  panel  method  can  calculate  the  vortex  bursting  of  the  vortex 
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iwake.  The  weighting  values  of  the  Kutta  equations  are  investigated 
for  the  leading-edges,  and  are  determined  to  be  problem  specific.  An 
^alternative  approach  to  the  management  of  the  wake  data  is  presented 
jand  implemented.  The  delta  wing  is  then  analyzed  for  20.5  and  30 
.degrees  angle  of  attack.  The  panel  method  produced  a  perturbation  and 
jenlargement  in  the  wake  at  the  two-thirds  chord  position  over  the  wing 
lat  30  degrees  angle  of  attack,  indicative  of  vortex  breakdown.  How- 
jever,  a  qualitative  comparison  with  the  results  of  the  20.5  degree 
jangle  of  attack  calculation  negates  this  conclusion.  Vortex  bursting 
.not  evident — over  the  wing  at — 30  degrees  angle  of  . - 

iVortex  Bursting,  Vortex  Breakdown,  ^6 — ^ — — - 

iVortex  Panel  Methods,  Vortex  Lattice  Methods 
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y(x)  -  ±^^fo.2969Vx-0.1260x-0.3516x^ 

0.2  V  (2-41) 

+  0. 2843x^-0. 1015x^ 

where , 

X  -  0.5 -^O.5cos(0)  (2-42) 

Where  0  is  incremented  from  0  to  2n  by: 

2ll  (2-43) 

A0  »  - - - 

number  of  panels 

This  was  used  to  capture  the  increased  curvature  of  the 
leading  edge. 

The  airfoil  was  modeled  using  200  panels,  at  an  angle 
of  attack  of  5  degrees,  and  was  impulsively  started.  The 

chord  was  set  equal  to  one,  and  the  time  step  was  0.01, 
since  this  was  the  approximate  average  length  of  the  panels. 
Figures  5  and  6  compare  the  vortex  sheet  at  time  steps  of 
25,  50,  75,  and  100.  The  vortex  cores  are  not  shown. 

Instead  a  curve  connecting  consecutive  cores  is  provided  as 
a  representation  of  the  vortex  sheet.  Notice  there  is  no 
visible  difference  in  the  wake  for  these  times.  Likewise, 
Figures  7  and  8  show  the  position  of  the  wake  at  125  and  150 
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time  steps.  A  small  difference  is  noticeable  at  150  time 
steps.  Both  algorithms  produced  essentially  the  same 
results.  However,  the  enhanced  splitting  method  required 
approximately  25  percent  more  computational  time.  This 
added  computational  requirement  would  increase  as  the 
calculation  was  continued.  The  primary  cause  was  the  added 
calculations  required  to  compute  the  local  velocity  twice 
for  the  entire  wake  each  time  step.  The  enhanced  splitting 
algorithm  therefore  did  not  provide  a  substantial  improve¬ 
ment  to  warrant  its  use  in  the  three-dimension  paneling 
method . 


Figure  5.  Wake  Formation  Using  the  Basic  Splitting  Algorithm 
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TS=25 


TS=50 


TS=75  TS=100 


NACA0012  Airfoil,  200  Panels,  5  Degree  Angle  of  Attack 

Figure  6.  Wake  Formation  Using  the  Enhanced  Splitting 
Algorithm 


7.  Downstream  Wake  Formation  Using  the  Basic 
Splitting  Algorithm 


2.00  2.10  '  2.20  . 2.30  2.40 . l.ho 

Enhanced  Splitting  Method 

NACA0012  Airfoil,  200  Panels,  5  Degree  Angle  of  Attock 


Figure  8.  Downstream  Wake  Formation  Using  the  Enhanced 
Splitting  Algorithm 


chapter  III 


Three-dimensional  Model 


3.1  Introduction 

A  three-dimensional  vortex  panel  method  was  used  to 
investigate  the  vortex  generated  over  a  unit  area  aspect 
ratio  delta  wing  at  high  angle  of  attack.  The  method  was 
developed  by  Mracek  and  uses  triangular  elements,  or  panels, 
on  which  the  surface  vorticity  is  piecewise  linearly  vary¬ 
ing.  Quadratical 1 y  varying  vortex  cores  are  placed  at  the 
edges  in  order  to  generate  a  nonzero  vorticity  at  the  edges. 
The  wake  was  modeled  as  a  vortex  sheet,  composed  of  discrete 
vortex  cores,  which  emanate  from  the  edges  of  the  surface. 
(5:Chapter  2) 

The  presentation  of  the  complete  aerodynamic  model  is 
too  extensive  to  present  here.  The  solution  method  used  was 
a  minimization  in  the  square  of  the  errors.  This  permited 
the  use  of  weighting  values  to  accentuate,  or  deemphasize 
certain  equations  with  respect  to  the  other  equations.  The 
weighting  value  for  the  Kutta  equations  was  investigated  for 
the  leading  edges. 

An  alternate  approach  to  Mracek's  modeling  of  the  wake 
was  performed.  The  new  method  is  analogous  to  a  finite 
element  technique,  using  nodes  and  elements.  The  nodes  are 
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convected  at  the  local  particle  velocity.  The  splitting 
algorithm  was  extended  from  the  two-dimensional  model  and 
applied  to  the  three-dimensional  wake. 

3.2  Aerodynamic  Model 

The  three-dimensional  vortex  panel  method  was  provided 
by  Mracek.  The  method  uses  triangular  panels,  or  elements, 
on  which  the  surface  vorticity  is  piecewise  linearly  vary¬ 
ing,  to  model  the  lifting  surface.  The  vertex  of  each 
triangular  element  was  termed  a  node.  Additional  vortex 
cores  are  placed  at  the  edges  in  order  to  generate  a  nonzero 
vorticity  at  the  edges.  The  wake  was  modeled  as  a  vortex 
sheet,  composed  of  discrete  vortex  cores,  or  vortex  lattice, 
which  emanated  from  the  edges  of  the  wing. 

Typical  to  vortex  paneling  methods,  there  exists  more 
equations  than  unknowns.  The  unknowns  to  the  problem  are 
the  strength  of  the  surface  vorticity  at  the  nodes,  and  the 
constants  of  integration  associated  with  the  surface  edge 
cores.  Again,  Laplaces  equation  was  satisfied,  and  only  the 
boundary  conditions  were  addressed.  The  resulting  matrix 
equation  was: 
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(3-1) 
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This  equation  is  the  culmination  of  four  sets  of  equations, 
each  of  which  are  discused  separately. 

The  no  penetration  conditions  at  the  midpoint  of  each 
element,  or  control  points,  is  given  by: 

V-n  *  O  on  the  surface  (3-2) 

Each  element  has  an  influence  in  the  velocity  at  every 
control  point  on  the  surface.  The  result  was  an  influence 
matrix,  [A]. 


/n\  (3-3) 

[A-C  D]|^|  -  <U-W> 

Also  included  was  the  influence  of  the  edge  cores  which  are 
dependent  on  the  surface  vorticity,  and  the  velocity  induced 
by  the  wake. 

The  flow  field  must  also  satisfy  the  continuity  equation 
for  incompressible  flow. 
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div  V 


O 


(3-4) 


A  vector  field.  A,  is  chosen  such  that 

V  *  curl  A 

Thus,  the  continuity  equation  becomes: 

div  (curl  A)  =  0 

The  vorticity  in  the  flow  field  is  defined  as: 

D  *  curl  V 

Because  the  freestream  is  i rrotational ,  the  vorticity 
defined  as  the  curl  of  the  disturbance  velocity  only, 
disturbance  velocity  can  then  be  found  by  solving  the 
following  equation: 


n  -  -v^A 


Green’s  theorem  yields: 


A(r.t) 


4  It 


!  r  -  s: 


(3-5) 


(3-6) 


(3-7) 
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The 


(3-8) 


(3-9) 
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From  Karamcheti,  if  the  vorticity  is  contained  within  a 
short  distance,  e,  of  the  surface,  as  e  approaches  zero,  the 
product  of  ln|  and  e  remains  constant,  and  D.  becomes  tangent 
to  the  surface  (3:530-532).  Therefore, 


•Y(s)=  _lim  He  (3-10) 

Rewriting  the  volume  integral  as  a  surface  integral: 


A(F) 


4nJ  ./s  I  r  -  s 


(3-11) 


Since  the  disturbance  velocity  is  the  curl  of  the  vector 
field  A,  and  the  divergence  of  the  curl  of  a  vector  is 
identically  zero,  any  vector  field  that  represents  the 
surface  vorticity  will  satisfy  the  continuity  equation  as 
long  as: 


divy  *  0  (3-12) 

This  divergenceless  condition  for  the  surface  can  be 
represented  in  matrix  form  by: 


[B 


{0> 


(3-13) 
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Quadratical 1 y  varying  cores  are  added  to  the  edges  of 
the  lifting  surface  to  either  capture  or  create  the  vortex 
sheet.  The  circulation  around  the  core  is  dependent  on  tne 
vorticity  of  the  neighboring  vortex  element.  This  is 
depicted  in  Figure  9. 


Figure  9.  Edge  Core  and  Edge  Element 
Along  the  edge  of  the  element: 


Yy2-  Yyl 

Y(X)  -  -X'Yy. 


(3-14) 


For  the  edge  cores  to  be  compatable  with  the  vortex  panel, 
the  circulation  is  given  by: 
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dr 

dx 


Y 


(3-15) 


Integration  provides  the  circulation  around  the  core. 


r(x) 


Vy2-Vyl 

d 


Yy,x  +  G 


(3-16) 


where  G  is  a  constant  of  integration.  In  order  to  make  the 
circulation  continuous  in  strength  around  the  body,  two 
adjoining  edge  cores  must  have  the  same  circulation  at  the 
connection.  Additionally,  the  circulation  around  the  cores 
are  continuous  through  the  first  derivative  at  the  juncture. 
These  compatibility  equations,  in  matrix  form,  are  given  by: 
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(3-17) 


The  final  matrix  equation  is  the  Kutta  equation.  The 
Kutta  condition  was  imposed  only  at  the  nodes  of  the  surface 
along  the  edges.  This  was  given  by: 


[K 


<P> 


(3-18) 
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A  more  in-depth  derivation  is  provided  in  the  following 
section . 

A  weighting  matrix  was  used  to  emphasize  one  equation 
with  respect  to  the  other  equations.  Hence, 
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(3-19) 


The  weighting  value  applied  to  the  divergenceless 
condition  has  little  impact  on  the  solution.  The  weighting 
value  used  for  this  condition  was  50.  The  weighting  value 
for  the  edge  compatibility  needed  to  be  large  because  these 
conditions  need  to  be  imposed  exactly.  Therefore,  the 
weighting  value  associated  with  the  edge  compatibility  was 
set  to  500.  Along  the  trailing  edge,  the  Kutta  condition 
also  needed  to  be  satisfied.  The  weighting  valued  applied 
to  the  Kutta  equation  was  set  to  50.  An  investigation  into 
the  weighting  value  of  the  Kutta  condition  along  the  leading 
edge  is  presented  in  section  3.3. 

Equation  3-19  can  be  expressed  as: 
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(3-20) 
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where , 
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The  chosen  method  of  solution  was  to  minimize  the  sum 
of  the  squares  of  the  errors.  This  is  the  solution  of  the 
following  equation: 


-  [S]^<T} 


(3-23) 


3.3  Kutta  Weighting  Values 

The  Kutta  condition  states:  (6:45-50) 


,,,,  -  O 


(3-24) 
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where  the  pressure  coefficient  is  defined  as: 


C, 


P-P. 


(3-25) 


Using  Bernoulli's  equation  for  unsteady  flow: 
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H(t) 


(3-26) 


where  R  is  the  position  vector  of  the  point  in  an  inertial 
reference  frame. 

The  far-fieid  boundary  condition  was  that  the  velocity 
potential  was  not  changing  with  respect  to  time.  The 
pressure  coefficient  can  then  be  written  as: 
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The  Kutta  equation  was  imposed  ideally  at  the  trailing 
edge.  This  was  performed  by  using  a  weighting  value  of  50 
for  the  Kutta  equations  associated  with  the  trailing  edge 
nodes.  Along  the  leading  edges,  however,  a  large  pressure 
gradient  exists.  In  order  to  capture  this  pressure  gradient 
a  large  number  of  elements  would  be  required.  This  was 
unachievable  along  the  front  of  the  wing  using  discrete 
points.  Therefore  the  Kutta  condition  can  not  be  satisfied 
ideally.  Instead,  since  it  was  known  that  the  wake  will 
convect  off  the  leading  edge  tangential  to  the  wing  surface, 
the  Kutta  condition  was  only  imposed,  in  as  much,  to  cause 
this  effect.  A  weighting  value  was  implemented  on  the  Kutta 
equations  associated  with  the  leading  edge  nodes. 

It  was  evident  that  the  required  weighting  value  for 
the  Kutta  equations  was  job  dependent.  Several  unit  area 
aspect  ratio  delta  wings,  modeled  with  various  number  of 
elements  were  analyzed.  Initially,  a  delta  wing,  modeled 
with  20  rows  of  elements  was  investigated  with  Kutta 
weighting  values  of  0.2  and  0.5.  The  angle  of  attack  was 
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20.5  degrees,  and  the  computations  were  computed  to  20  time 
steps,  or  one  chord  length  of  wake.  The  Kutta  weighting 
value  of  0.2  was  not  sufficient  enough  to  establish  the 
tangency  condition.  The  value  of  0.5  however,  produced  the 
satisfactory  result.  As  a  test  to  investigate  if  the  0.5 
value  was  too  excessive,  another  computation  was  performed 
using  0.45.  This  value  proved  inadequate,  therefore,  a 
calculation  with  a  Kutta  weighting  value  of  0.48  was 
performed  producing  the  required  tangency  condition.  A  unit 
area  aspect  ratio  delta  wing  was  then  modeled  with  30  rows 
of  elements,  using  the  0.48  Kutta  weighting  value.  Again 
the  problem  was  carried  out  to  obtain  one  chord  length  of 
wake,  or  in  this  case,  30  time  steps.  For  this  calculation, 
the  value  of  0.48  was  inadequate.  The  required  Kutta 
weighting  value  for  the  30  row  delta  wing  was  determined  to 
be  0.5.  The  angle  of  attack  was  then  increased  to  30 
degrees.  Initially  a  value  of  0.5  was  used  as  the  weighting 
value.  This  proved  inappropriate  and  a  value  of  0.515  was 
needed  for  this  increased  angle  of  attack.  No  analytical, 
or  empirical  relationship  could  be  derived  to  determine  the 
appropriate  weighting  value  for  any  job  in  general.  This 
data  is  provided  in  Table  1. 
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Table  1. 


Kutta  Weighting  Value 


No. 

Angle  of 

Kutta 

Resul ts 

of 

ft  1 1  a  C  V 

W(»:i  "i  ng 

Rows 

Value 

20 

20.5 

0.200 

Insufficient 

20 

20 . 5 

0 . 500 

Over -cons trained 

20 

20 . 5 

0.450 

Insufficient 

20 

20.5 

0.480 

Satisfactory 

30 

20 . 5 

0.480 

Insufficient 

30 

20.5 

0.500 

Satisfactory 

30 

30.0 

0.500 

Insufficient 

30 

30.0 

0.515 

Satisfactory 

3.4  Wake  Modeling 

The  wake  was  model  as  a  vortex  sheet  comprised  of 
discrete  vortex  cores.  The  fluid  flow  was  started  impul¬ 
sively,  therefore,  initially  there  was  no  wake.  Once  the 
fluid  was  set  into  motion,  vorticity  was  formed  on  the 
surface.  With  this  quadratical 1 y  varying  vortex  cores  were 
developed  along  the  edges.  The  movement  of  the  wake  was 
given  by  the  no  pressure  jump  condition.  The  pressure  on 
the  upper  surface  of  the  wake  must  equal  the  pressure  on  the 
lower  surface  of  the  wake.  In  order  to  satisfy  this  condi¬ 
tion,  the  wake  cores  were  convected  at  the  local  particle 
velocity.  The  cores  new  positions  were  given  by: 


rCt  +  At)  -  r(1:)-<-VpAt 


(3-29) 
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where , 


Vp  =  local  particle  velocity 
dt  =  time  step 

Once  the  cores  have  moved  downstream,  the  flow  field  around 
the  wing  is  changed.  This  required  a  new  vortex  distribu¬ 
tion  on  the  wing  and  edges.  The  strength  of  the  wake  was 
determine  by  the  conservation  of  circulation. 


This  procedure  was  performed  throughout  the  calculation  for 
each  time  step. 

In  order  to  simplify  the  wake,  vortex  cores  are  ai  o 
used  to  model  the  wake  parallel  to  the  flow.  The  endpoints 
of  which  connect  to  the  endpoints  of  the  cores  convected 
from  the  edges  of  the  surface.  This  creates  a  quadrilat¬ 
eral,  or  wake  element,  consisting  of  four  vortex  cores 
connected  at  their  endpoints.  Figure  10,  illustrates  the 
development  of  the  wake  lattice.  Where  the  cores  connect 
are  termed  wake  nodes.  To  satisfy  the  fact  that  the  vortex 
cores  that  comprise  the  wake  will  be  of  constant  strength, 
the  core  circulation  will  be  the  average  value  of  the 
guadratical ly  varying  edge  core.  The  quadrilateral  bounds  a 
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closed  ring  of  circulation.  Only  the  value  of  this 
circulation  need  be  stored.  The  circulation  of  the 
quadratical ly  varying  edge  core  is  given  by: 


r(x) 


o 

+  g3 


(3-31) 


Therefore,  the  strength  of  the  ring  of  circulation  is  given 
by: 


r(x)dx 

Jo _ 

ci 


(3-32) 


As  the  calculation  is  carried  out,  the  wake  moves 
downstream  and  spreads  out.  As  with  the  two-dimensional 
calculations,  the  wake  was  split  in  order  to  more  closely 
model  a  continuous  vortex  sheet.  The  model  was  an  extension 
of  the  one  chosen  in  the  two-dimensional  model,  although, 
now  the  wake  was  split  in  the  spanwise  direction.  Once  a 
wake  element  has  deformed  to  warrant  splitting,  it  was 
divided  into  two  new  wake  elements,  and  two  wake  nodes  are 
added.  This  is  represented  in  Figure  11.  The  strength  of 
the  new  wake  elements  are  given  by: 
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Figure  11.  Splitting  the  Wake  Lattice 


Mracek  performed  the  same  splitting  technique,  but 
managed  the  wake  data  in  a  redundant  manner.  He  considered 
each  row  of  wake  data  as  an  independent  wake  entity,  which 
necessitated  the  existence  of  two  rows  of  wake  nodes  for 
each  row  of  wake  elements.  To  make  the  wake  continuous,  the 
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rows  of  wake  elements  would  be  attached.  Assuming  there 
were  the  same  number  of  wake  elements  in  two  consecutive 
rows,  the  upstream  row  of  wake  nodes  would  correspond  to  the 
downstream  row  of  wake  nodes  from  the  preceding  row  of 
elements.  This  redundancy  required  twice  the  amount  of 
computer  memory  needed  to  describe  the  wake. 

To  eliminate  the  extra  computer  memory  requirement,  the 
bookkeeping  of  the  wake  used  a  finite  element  technique. 

The  poLtions  of  the  nodes,  and  the  strengths  of  the 
circulations  of  the  wake  elements  were  only  stored.  Four 
nodes  are  associated  with  each  wake  element.  Each  node  is 
shared  by  the  surrounding  elements.  The  wake  can  then  be 
represented  two-dimensional ly.  The  nodes  are  numbered 
sequentially,  starting  at  the  upper  left  side  and  ending  at 
the  bottom  'ight  side  of  the  wake.  A  matrix  consisting  of 
rows  and  columns  was  used  to  manage  the  circulation  strength 
of  the  wake  elements.  Figure  12. 

A  two-dimensional  representation  of  the  wake  is  given 
by  Figure  12.  For  the  next  time  step  this  wake  was  used  to 
in  the  solution  process  to  solve  for  the  vorticity 
distribution  over  the  wing.  It  was  also  used,  along  with 
the  new  vorticity  distribution,  to  determine  the  local 
particle  velocity  of  the  wake  nodes.  The  wake  nodes  were 
then  convected  at  the  local  particle  velocity.  A  new  row  of 
wake  nodes  was  added  along  the  edges  of  the  wing  along  with 
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Figure  12.  Finite  Element  Representation  of  the  Wake 


a  new  row  of  wake  elements.  The  spanwise  distance  of  the 
elements  was  then  checked  to  see  if  splitting  was  necessary. 
Additional  nodes  and  elements  were  added  as  required.  This 
process  starts  with  the  last  row  of  elements  and  proceeded 
upstream  towards  the  wing.  Once  a  row  of  elements  was 
reached  that  does  not  require  splitting  and  still  contained 
its  original  number  of  elements,  no  further  investigation 
into  the  spanwise  distance  of  the  elements  was  needed.  This 
was  based  on  the  premise  that  the  wake  is  continually 
diverging.  The  nodes  and  elements  were  then  resequenced  to 
conform  to  the  previous  numbering  scheme,  Figure  13. 

As  a  check,  two  calculations  using  a  unit  aspect  ratio 
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delta  wing,  modeled  with  20  rows  of  elements,  at  an  angle  of 
attack  of  20  degrees  were  performed  each  using  one  of  the 
wake  data  management  systems  described  above.  The  calcula¬ 
tions  were  performed  to  20  time  steps.  The  finite  element 
methodology  was  verified,  computing  the  same  wake  mesh  as 
Mraceks  method.  However,  the  finite  element  technique 
required  less  total  memory,  and  computer  time  than  the 
previous  method  employed.  Therefore,  this  wake  data 

management  procedure  was  retained  for  the  remainder  of  this 
effort . 


Figure  13.  Three-Dimensional  Wake  Formation 


chapter  IV 


Resul ts 


Two  geometries  were  investigated  for  the  analysis.  In 
both  cases  a  unit  area  aspect  ratio  delta  wing  was  modeled 
using  30  rows  of  elements.  The  wing  was  analyzed  at  20.5 
and  30  degrees  angle  of  attack.  It  is  known  that  vortex 
breakdown  does  not  occur  over  a  delta  wing  with  an  aspect 
ratio  of  one  at  20.5  degrees  angle  of  attack.  Therefore, 
this  was  calculated  as  a  baseline.  An  angle  of  attack  of  30 
degrees  was  used  to  determine  if  paneling  methods  can  com¬ 
pute  the  vortex  bursting  that  occurs  over  the  delta  wing. 

In  order  to  establish  if  vortex  bursting  occurs,  the 
criteria  established  by  O'Neil,  Barnett,  and  Louie  shall  be 
used.  In  both  cases,  the  fluid  was  impulsively  started. 

The  Kutta  weighting  values  for  the  20.5  and  30  degree  angle 
of  attacks  were  0.50  and  0.515  respectively.  These  values 
were  sufficient  to  sustain  the  wake  tangency  condition  with 
the  surface  of  the  airfoil  throughout  the  calculations. 
Figure  14  presents  the  element  composition  of  the  delta 
wing . 

Both  calculations  were  carried  out  to  produce  one  chord 
length  of  wake.  The  time  step  was  set  equal  to  the  height 
of  the  triangular  elements.  This  produced  wake  elements 
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Figure  14.  Element  Composition  of  a  Unit  Area  Aspect  Ratio 
Delta  Wing 


that  were  initially  geometrically  square  and  required  30 
time  steps  to  produce  one  chord  length  of  wake.  Figures  15, 
and  16  show  the  wake  mesh  at  this  time  for  the  20.5,  and  30 
degrees  angle  of  attack.  The  view  was  from  above  the  air¬ 
foil.  In  both  figures  the  wing  elements  were  not  presented 
to  help  clarify  the  pictures.  The  'X'  represents  the  body 
frame  origin,  and  the  front  apex  of  the  wing.  The  trailing 
edge  is  easily  recognized  as  the  vertical  line  where  the 
wake  is  orginating.  The  paneling  method  computes  both  the 
leading  edge  and  the  trailing  edge  vortices  typical  of  the 
flow  over  delta  wings.  However,  the  method  breaks  down  for 
the  30  degree  angle  of  attack  where  these  two  vortices  meet 
and  combine.  This  was  due  to  the  first  order  convection 
scheme  used  for  the  wake  nodes,  convecting  the  cores  at  the 
local  particl  velocity.  As  two  vortex  cores  come  within 
close  proximity  of  each  other,  their  respective  velocities 
increase  dramatically.  This  is  evident  from  equation  2-29. 
Figures  17,  and  18,  provide  a  hidden  line  picture  that  pro¬ 
vides  a  clear  view  of  the  leading  edge  vortices  and  the 
destruction  of  the  wake  mesh  for  the  30  degree  angle  of 
attack.  For  the  20  degree  angle  of  attack,  the  leading  edge 
vortex  expands  in  a  continuous  manner  as  it  is  convected 
downstream.  However,  viewing  the  leading  edge  vortex  for 
the  30  degree  angle  of  attack,  shows  a  perturbation  and 
enlargement  of  the  vortex  approximately  at  the  two-thirds 
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Figure  15.  Wake  Mesh  at  20.5  Degrees  Angle  of  Attack  and  30 
Time  Steps 


Figure  16.  Wake  Mesh  at  30  Degrees  Angle  of  Attack  and  30 
Time  Steps 
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cord  position.  This  fits  the  description  of  vortex- 
breakdown  provided  by  O’Neil,  Barnett,  and  Louie.  Cross- 
sections  of  the  leading-edge  vortices  are  taken  in  the  Y-Z 
plane  just  along  the  trailing  edge  of  the  wing.  Figures  19, 
and  20.  Here  the  paneling  method  performed  well  by  modeling 
the  vortex  roll-up.  The  vortex  for  the  30  degree  angle  of 
attack  is  larger  than  the  vortex  produced  at  an  angle  of 
attack  of  20  degrees,  as  would  be  expected.  Further  inves¬ 
tigation  into  the  pressure  gradients  in  the  spanwise  direc¬ 
tion  should  be  performed  to  verify  this  result. 

Figures  21,  and  22,  show  the  coefficient  of  pressure 
distribution  in  the  spanwise  directions  for  20.5  and  30 
degrees  angle  of  attack  respectively.  The  length  in  the 
spanwise  direction  is  normalized  by  one-half  the  local  span 
(Ymax) .  The  pressure  distributions  are  provided  at  one- 
half,  two-thirds,  and  five-sixths  chords  and  at  the  trailing 
edge.  These  curves  are  typical  of  the  trends  investigated 
experimentally  under  similar  conditions  over  delta  wings. 

It  is  evident  here  that  the  Kutta  condition  is  not  satisfied 
at  the  leading  edges  as  the  coefficient  of  pressure  does  not 
equal  zero  at  the  end  of  each  curve,  as  it  should.  The  rise 
in  the  curves  are  the  enhanced  lift  provided  by  the  suction 
induced  by  the  leading  edge  vortices.  A  peculiar  decrease 
in  the  coefficient  of  pressure,  or  decrease  in  lift,  is  con¬ 
spicuous  at  the  two-thirds  chords  position  at  30  degrees 
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Figure  19.  Trailing  Edge  Cross  Section  of  the  Wake  for 
20.5  Degrees  Angle  of  Attack 


Figure  20.  Trailing  Edge  Cross  Section  of  the  Wake  for 
30  Degrees  Angle  of  Attack 
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Pressure  Coefficient  Pressure  Coefficient 
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Figure  21.  Pressure  Distribution  Over  the  Wing  (Cp) 
Angle  of  Attack  =  20.5  degrees 
One  Chord  Length  of  Wake 


angle  of  attack.  This  corresponds  to  the  chordwise  location 
of  the  vortex  perturbation.  These  sharp  gradients  in  the 
pressure  coefficient,  however,  contradict  the  vortex  break¬ 
down  criteria  as  stated  previously. 

In  an  attempt  to  eliminate  the  sharp  pressure 
gradients,  the  calculations  were  continued  for  another  ten 
time  steps,  or  one-third  chord  length.  This  was  performed 
to  allow  the  wake  to  split  more  thoroughly,  and  hence, 
create  the  mild  pressure  gradients  required  to  meet  the  vor¬ 
tex  breakdown  criteria.  Figures  23,  and  24,  show  the  wake 
mesh,  at  this  time,  for  the  20.5  and  30  degrees  angle  of 
attack  respectively.  At  this  time  the  wake  is  beginning  to 
deteriorate  for  the  20.5  degrees  angle  of  attack.  The  wake 
has  severely  degenerated  for  the  30  degree  angle  of  attack 
behind  the  wing,  although,  it  still  retains  its  structure 
over  the  wing.  The  perturbation  in  the  leading-edge  vortex 
is  still  evident  and  implies  the  existence  of  vortex  break¬ 
down  , 

Again  the  pressure  distributions  in  the  spanwise 
direction,  at  the  same  chordwise  positions,  over  the  wing 
are  reviewed.  The  coefficient  of  pressure  distributions  are 
provided  in  Figures  25,  and  26.  Each  position  along  the 
surface  of  wing  at  20.5  degrees  angle  of  attack  provides 
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Figure  23.  Wake  Mesh  at  20.5  Degrees  Angle  of  Attack  and 
40  Time  Steps 


X 


Figure  24.  Wake  Mesh  at  30  Degrees  Angle  of  Attack  and 
40  Time  Steps 
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more  lift  with  the  increased  wake.  The  pressure  gradient  at 
the  two-thirds  chord  position,  at  an  angle  of  attack  of  30 
degrees,  has  not  diminished,  and  has  actually  increased  with 
the  increased  amount  of  wake.  However,  this  is  the  only 
position  that  exhibits  a  decrease  in  the  amount  of  lift  pro¬ 
vided  for  both  calculations.  This  severe  pressure  gradient 
again  fails  the  criteria  for  vortex  bursting  over  the  wing. 

Due  to  the  extreme  wake  degeneration  for  the  30  degrees 
angle  of  attack,  a  question  arises  if  this  was  a  valid  cal¬ 
culation.  In  an  attempt  to  answer  this  query,  another  cal¬ 
culation  was  performed  for  the  30  row  delta  wing  at  an  angle 
of  attack  of  30  degrees.  However,  for  this  computation  the 
wake  was  not  allowed  to  split  beyond  the  one  and  one-third 
chords  position  from  the  leading  apex  of  the  wing.  The  wake 
mesh  is  provided  in  Figure  27.  For  this  calculation  the 
wake  has  not  deteriorated  to  the  point  of  the  previous  com¬ 
putation.  However,  the  coefficient  of  pressure  distribution 
over  the  wing  is  equivalent  to  the  preceding  coefficient  of 
pressure  distribution.  Figure  28.  This  provides  the  unex¬ 
pected  result  that  the  pressure  distribution  over  the  wing 
is  relatively  independent  to  the  wake  downstream. 

A  comparision  of  the  wake  cross  sections  at  two-thirds, 
and  five-sixths  chords  length  and  at  the  trailing  edges  are 
provided  for  the  30  and  20.5  degrees  angle  of  attack 
computations  in  Figures  29,  30,  and  31,  respectively.  The 
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Figure  27.  WaKe  Mesh  at  30  Degrees  Angle  of  Attacjc  and 
40  Time  Steps,  No  Splitting  Beyond  -40 
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comparisons  are  performed  at  40  time  steps.  At  the  two- 
thirds  chord  position,  the  wake  has  only  progressed  one  com¬ 
plete  revolution.  At  the  trailing  edge,  it  has  developed  to 
approximately  one  and  one-half  revolutions.  At  every 
chordwise  postion,  the  wakes  generated  by  the  two  calcula¬ 
tions  are  analogous.  The  prominent  distinction  between  the 
different  angle  of  attacks  is  the  magnitude  of  the  leading 
edge  vortex.  Again  this  leads  to  the  conclusion  that  vortex 
breakdown  does  not  exist  in  the  30  degrees  angle  of  attack 
calculation. 
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Figure  29.  Vortex  Comparisons  At  the  Two-Thirds 
Chordwise  Position 
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Chapter  V 


Conclusions  and  Recoiromendations 


A  unit  aspect  ratio  delta  wing  was  analyzed  to  deter¬ 
mine  if  a  unsteady  vortex  panel  method  could  compute  the 
bursting  of  the  leading  edge  vortex.  The  criteria  for  which 
vortex  bursting  would  be  identified  was  that  defined  by  O'N¬ 
eil,  Barnett,  and  Louie;  a  degeneration  of  the  leading-edge 
vortex  combined  with  relatively  mild  gradients.  (6:220)  A 
two-dimensional  vortex  panel  method  was  employed  to  evaluate 
wake  splitting  procedures  of  the  vortex  cores  representing 
the  vortex  sheet,  or  wake.  The  worthier  splitting  scheme 
was  then  extended  to  the  three-dimensional  panel  methods. 

The  weighting  values  applied  to  the  Kutta  equations  along 
the  leading  edges  were  determined  to  be  problem  dependent. 

No  relationship  between  the  airfoil  parameters  and  the 
weighting  value  could  be  determined.  A  finite  element  tech¬ 
nique  was  employed  to  manage  the  wake  data.  This  provided  a 
reduction  in  the  amount  of  computer  time  and  memory  required 
for  a  given  calculation.  A  unit  area  aspect  ratio  delta 
wing  analysis  was  performed.  A  qualitative  comparison 
between  20.5  and  30  degrees  angle  of  attack  was  accom¬ 
plished.  At  the  30  degrees  angle  of  attack  the  paneling 
method  produced  a  perturbation  and  enlargement  in  the 
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leading  edge  vortex  at  the  two-thirds  r-hordwise  position. 
However,  wake  cross-sections  and  coefficient  of  pressure 
comparisons  to  the  20.5  degrees  angle  of  attack  computation 
determined  that  the  panel  methods  did  not  calculate  vortex 
breakdown . 

A  possible  explanation  of  the  failure  to  predict  vortex 
bursting  was  the  splitting  method  employed.  The  process 
would  check  the  spanwise  distance  of  the  wake  elements,  and 
split  the  elements  in  the  spanwise  direction  if  warranted. 
This  was  based  on  the  premise  that  the  wake  is  continuously 
diverging.  However,  the  wake  elements  would  also  elongate 
in  the  chordwise  direction.  This  is  apparent  in  Figures  29 
through  31.  The  area  most  susceptible  to  this  stretching  is 
in  the  leading  edge  vortex.  A  method  of  splitting  the  wake 
in  the  chordwise  direction  should  be  performed  before  fur¬ 
ther  analysis  of  vortex  breakdown  is  accomplished. 

The  vortex  sheet  generated  off  the  edges  of  the 
airfoil,  in  this  effort,  was  used  as  the  flow  visualization 
of  the  wake.  Further  analysis  could  be  performed  by  the 
inclusion  of  Lagragian  points  within  the  flow  field.  This 
could  easily  be  incorporated  within  the  computer  code  and 
would  possibly  provide  further  insight  into  the  problem. 

Another  area  of  interest  would  be  the  implementation  of 
a  second  order  convection  scheme  of  the  wake  nodes.  The 
intermingling  of  the  leading  edge  and  trailing  edge  vorti- 
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cies  induced  unrealistically  high  velocities  of  the  wake 
nodes.  Thus,  the  deterioration  of  the  wake  mesh.  A  second 
order  convection  scheme  could  possibly  eliminate  this  defi¬ 
ciency.  However,  the  use  of  a  higher-order  method  would 
require  more  computer  time,  but  should  still  be  performed  to 
determine  its  benefits. 
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